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We report non-resonant x-ray Raman scattering (XRS) measurements from hexagonal boron 
nitride for transferred momentum from 2 to 9 A along directions both in and out of the basal 
plane. A symmetry-based argument, together with real-space full multiple scattering calculations 
of the projected density of states in the spherical harmonics basis, reveals that a strong pre-edge 
feature is a dominantly Yio-type Frenkel exciton with no other s-, p-, or d- components. This 
conclusion is supported by a second, independent calculation of the q-dependent XRS cross-section 
based on the Bethe-Salpeter equation. 
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The physics of low-energy photoelectrons in solids is 
a complex, many body problerrfl'. All aspects of the 
electro nic structure of the material must be taken into 
accounlP^, including the interaction betwee n the photo- 
electron and the complementary hol^^EIIIll This latter 
effect allows for similar long-lived bound states in metals 
(generally called Fermi-edge singularities)^ and in insula- 
tors (called excitons)!^. As a canonical example of a many 
body problem in condensed matt er ph ysics, core-excitons 
are a topic of continuing interestP'^. 

The present synergy betw een stea dily progressing ah 
initio theoretical treatmentd^l^lllini and on-going im- 
provements in instrumentation^^ in studies of non- 
resonant x-ray Raman scattering (XRS) shows a s trong 
potential for rapid progress on this old problerrP^. XRS 
is the inelastic scattering of hard x-rays from bound elec- 
trons, and XRS studies have seen a recent explosion in 
the number and range of applicationj^' ^ ' " ' ^^ ' ^^ ' ^^ ' ^^ -^Sim. 
In comparison with q-resolved electron energy loss spec- 
troscopy (EELS)liS!, XRS is more suited for bulk con- 
densed matter systemJ^^ due to the large penetration 
length afforded by the relative high-energy incident and 
scattered photons. The measured q-dependent XRS pro- 
vides direct information about multipol e contribut ions 
to the dynamic structure factor S{(i,uj^^^^. It 
can be d irectly co mpared with theoretical calculations 
of S (q, |^)| 6 | i4 | i5 [i6]^ g^jj^j jj-^ some cases be inverted to 
provide an experimental meas ure of the perturbed pro- 
jected density of states (Z-D0S)pE2l^ which can again be 
compared to theory. 

XRS has shown itself to be especially suitable for spec- 
troscopy of the angular characteristics of core excitons at 
the near edge region. For example, XRS studies^ ^ on 
LiF and icosahedral B4C convincingly demonstrated s-, 
and p-type excitons at the F and B K-edges, respectively, 
which can be attributed to atoms located at a center of 
inversion symmetry of the unit cell with parity a good 



quantum number for the final stated. Here, we demon- 
strate that for some systems one can learn not only the 
A/ selection rule for the exciton but also a full descrip- 
tion of its angular characteristics in terms of spherical 
harmonics, l^m- 

We present measurements of the q-dependent XRS of 
hexagonal boron nitride (h-BN), an anisotropic system 
with all atoms sitting at positions with reflection symme- 
try to the basal plane (Figure 1, inset )^'^. The anisotropy 
of the dipole-limit core excitation spectra for q |j c and 
q _L c for h -BN is well known from previous EEL^^, 
XANES22123|^ and XRS^"* studies at the g = limit. We 
extend on this prior work by measuring the q-dependent 
boron XRS for both q || c and q _L c out to q = 9A , 
a momentum transfer that is clearly beyond the dipole 
limit for the B Is initial state. By expanding the final 
state wavefunction onto the spherical harmonics basis, 
we identify a predominantly yio-type exciton in the pre- 
edge region. This conclusion is critically supported by 
two independent ah initio calculations, one of the per- 
turbed /-DOS and the other a direct calculation of the 
q-dependent XRS. 

Our sample is highly-oriented compression-annealed 
pyrolytic h-BN (Advanced Ceramics Research with di- 
mensions 14 X 5 X 0.9 mm'^. The crystallite size is of 
the order of 1000 A with a mosaic spread about 2 to 3 
degrees along the c-axi^221 The orientation within the 
basal plane is assumed to be random. All measurements 
were performed at the XOR/PNC-20-ID beamline of the 
Advanced Photon Source at Argonne National Labora- 
tory. An 11.4-cm-diameter spherically bent Si (5, 5, 5) 
crystal analyzer was mounted in the vertical scattering 
plane on the arm of a Huber diffractometer, forming a 
near-backscattering Bragg diffraction geometry with the 
Amptek Si-PIN photodiode detector just above the sam- 
ple position. The incident x-ray energy was scanned with 
the x-ray analyzer energy fixed at 9890 eV. Henceforth, 
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FIG. 1; (color online). (Inset): The crystalline structure of 
h-BN has 2-dimension hexagonal sheets stacked directly on 
top of each other, with the boron (black) and nitrogen (gray) 
always being nearest neighbors^" . (Main figure) : The normal- 
ized XRS spectra from h-BN for q || c (top panel) and q _L c 
(bottom panel). The dashed line on the x-axis indicates the 
integration boundary for the pre-edge peak intensity shown 
in Fig. 4. 



we refer to all energies in terms of energy loss with respect 
to the elastic scattering peak. Inelastic x-ray scattering 
spectra were measured at q ~ 2.08,3.83,5.01,6.10,7.50, 
and 9.01 A for directions parallel and perpendicular to 
the c-axis, with a resolution varying from 0.25 to 0.11 

A . Further details of our apparatus have been reported 
elsewhere--". All of the elastic peaks are fit well by a 
Gaussian function, defining the energy resolution to be 
between 1.32 eV and 1.66 eV FWHM with the broaden- 
ing a consequence of the different effective sample size as 
a function of scattering angle. 

The separation of the XRS signal from the larger 
Compton background is a k ey experimental difficulty in 
all XRS measurement^iS'^II The valence Compton back- 
grounds are approximated by fitting the pre-edge data to 
standard Compton forms and using the low-q limit XRS 
form as a referencd^. The data were numerically broad- 
ened, as necessary, to the maximum resolution of 1.66 eV 
FWHM for visual comparison. 

The resulting boron XRS spectra are presented in Fig- 
ure 1. The pre-edge peak position is 191.9 ± 0.1 eV for 
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FIG. 2: (color online). Ab initio FEFF8 calculation of the l- 
DOS onto different spherical harmonics with odd (top panel) 
and even (bottom panel) parity for an h-BN cluster of 329 
atoms. See the text for details. 



the c-axis spectra, and the main edge peak position is 
198.8±0.1 eV along both directions, con sistent wit h prior 
work using other experimental method^^HHIlIIlll, All of 
the spectra were normalized according to the integrated 
intensity in the 210 eV to 250 eV range, mimicking the 
/-sum rule, and naturally introducing a factoiJ^. The 
resulting XRS spectra are shown in Fig. 1 for both q || c 
(top panel) and q _L c (bottom panel). 

Previous experimental studies have suggested that the 
physic al ori gin of the pre-edge peak in h-BN is a core 
excitor|22129| Theoretical calculations have confirmed this 
notion by demonstrating the importance of the interac- 
tion between the core-hole and the photoelectron"^" . In 
the XRS spectra for q || c (Fig. 1, top panel) this exci- 
ton state exhibits a prominent monotonically decreasing 
peak intensity with increasing q. Another less obvious 
but equally significant feature is the total absence of the 
exciton peak in spectra measured within the basal plane 
for any q values. Although the shape and intensity of 
the main edge peaks along both directions also display 
subtle variations, we concentrate on the behavior of the 
exciton. 

We now present a general formalism for q-dependent 
XRS in order to understand the interplay between the 
final-state characteristics measured by core-exciton spec- 
troscopy and the local symmetry at the sites of interest. 
The dynamic structure factor S'(q, uj) is expressed in a 
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FIG. 4: (color online), (solid dots with error bars, red): In- 
tegrated exciton intensities from the experimental XRS data 
for energies below 194.5 eV. (open squares, black): Theoreti- 
cally calculated exciton intensity (Fig. 3) integrated over the 
same energy range and plotted with an overall adjustable scal- 
ing factor to the experimental data, (dashed curve, blue): A 
simple calculation of S{q)/q^ using boron Is core-state and a 
crude hydrogen-like final state with a Ztoron + 1 approx- 
imation for the core-hol d^ ^ ^^ l 



FIG. 3: (color online). First-principles calculation of S{q^,u}) 
over the experimentally measured q range along both (a) q _L 
c and (b) q || c directions. See the text for details. 



single-electron picture as 

5(q,c.) = J2 l(/|e'''""|0)p5(AiJ - M- (1) 
/ 

Following the short-range order theor j^^S" due to the local- 
ized core-hole effect and the general guideline of Doniach 
et al^, the angular parts of both the electron wave func- 
tions and transition matrix element are expanded in the 
spherical harmonics basis as 

^/(r) = ^i?L,(r)y,,,„,(0,^) (2) 

I'm' 

and 

e^q-r = inY,^'ji{qr)Yira{Og,^g)Y;^{e,^), (3) 

l,m 

with R{i^i{r) = (Yi'm'l/). The initial core state is as- 
sumed to be s-type, and the final state is usually a mix- 
ture of partial waves because of the ligand-field influence 
in a solid system. Substitution of these expansions into 
Eq. 1 separates the angular variables of the transferred 
momentum q from its magnitude in S'(q, lo) as 

5(q, c^) = ^ I ^ SfUq)yim{0g, ^g)\^6{AE~hLo), (4) 



where 

/•oo 

SfUq) = i'V4^ jl{qr)R{^{r)Ro{r)-r^dr (5) 
Jo 

is the magnitude of the Yim-projected momentum space 
convolution between the initial state |0) and each final 
state IfP^. This approach is very similar to the formal- 
ism of Soininen et al^^, but without the directional av- 
eraging (i.e., summation over m) employed in that work. 

For the first three orders of spherical harmonics expan- 
sion (/ up to 2), partial waves of loo, ^i,±i, ^20, ^md 1^2, ±2 
could all have non- vanishing contributions to ^(q, to) for 
q _L c. Hence, the non-existence of the exciton peak 
for S'(q, Lu) measured along q _L c for any magnitude of 
q indicates the absence of these partial waves. Also the 
y2.±i do not have contributions to the XRS spectra along 
either direction. This leaves Yio the only spherical har- 
monic component for Z < 3 to contribute to the q || c 
spectra. 

This identification of a Yig-type exciton is consistent 
with the fact that boron atoms sit only at sites with 
refiection symmetry about the a-b basal plane, which 
demands that each final state be an eigenstate of the 
corresponding refiection parity^^. For a predominantly 
Yio-type exciton, other possible higher-order Yim com- 
ponents must possess the same refiection parity. This 
symmetry argument puts stringent restrictions on other 
possible spherical-harmonic components to the final state 
wave function. Note that Yqo, yi,±i, Y20, and 1^2, ±2 all 
have even parity about the a-b plane, while only Yiq and 
Y2.±i have odd parity. Given the additional refiection 
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symmetry about the a-c plane for boron sites (Fig. 1, 
inset), F2,±i could be also excluded from the final state 
wave function, because of the different a-c plane reflec- 
tion parity in comparison with Yiq. Thus for / < 3, the 
exciton state has a purely Yiq type of wave function. This 
is our central result. 

Compared to previous anisotropy studies by XANES 
and EELS that only excluded Yi^±i contributions be- 
cause of the dipole selection rul^^U^ the present q- 
dependent XRS study provides a more complete picture 
of the exciton state characteristics as of dominantly p^- 
type (^lo), with no s-, Px,y-, or d-components. The above 
analysis, in which the q-dependent XRS spectra may be 
qualitatively interpreted in the context of the constraints 
imposed by the crystal symmetry, should be generically 
applicable to the analysis of the angular characteristics 
of excitons in other high-symmetry crystals. 

Seeking independent verification and further insight 
into the exciton 's q-dependence, we performed two in- 
dependent first principles theoretical calculations. First, 
we modified the FEFF^^H code to calculate the density of 
final states projected onto the individual spherical har- 
monics as {Yim\G{r,r,uj)\Yim), with the core-hole effect 
included (Figure 2). The calculation confirms our parity 
conservation argument by showing no contribution from 
even reflection parity Yjm's and also indicates the Z-DOS 
for other odd Yim's with / up to 3 are negligibly small. 
This includes Y^q, which is the only other possible con- 
tributor to the q II c spectra for / up to 3. Note that 
the single Yiq characteristic of the wave-function indi- 
cates the exciton state is localized at the scattering atom. 
Hence it could not be a superposition of 2p states from 
both the excited boron atom and its neighboring nitrogen 
atoms, as typically in either a band structure hybridiza- 
tion or the molecular orbital theory. Such a superposition 
of wave functions from different atomic sites would not 
generate a single spherical harmonic wave function. This 
is a significant refinement of the prior understanding of 
this exciton peak as a Is — tt* transitiorl^illlllll 

Second, we calculated the XRS spectra from the boron 
K-edge using a computational Bethe-Salpeter equation 
method for modeling core-excited stateJ^S while includ- 
ing calculation of the Hedin-Lundqvist's GW self-energy^ 
for both the energy dependent quasi-particle shifts and 
the quasi-particle lifetime broadening. Convergence was 
checked with respect to the numerical cutoff^SH xhe sim- 
ulated spectra along both axes are plotted in Figure 3, 
using the same normalization and resolution as the ex- 
perimental data. The simulations also verify that the 
general spectral profile far away from the near edge has 
little variation in the probed q range. 

We plot in Figure 4 the integrated exciton intensity for 
energies below 194.5 eV as a function of g for q || c for 
both experimental data (solid dots with error bar) and 
the K-edge calculation (open squares). The estimated 
relative standard uncertainty of the integrated intensity 
is 4% from a combination of statistical and estimated 
systematic errors applying to each independent measure- 
ment. Due to the weak g-dependence of the atomic 



backgrouncPfil, our method of normalizing the spectra 
may also introduce another 20% systematic error into the 
overall, smooth g-dependence of the integrated area un- 
der the pre-edge feature. Hence the agreement between 
theory and experiment is indeed good, but certainly not 
so exacting as is shown in Fig. 4. As the theory has no 
lattice motion (phonon) included, the Frank-Condon ef- 
fect from electron-phonon couplinJ^ISl! jg not critical and 
the momentum is largely conserved during scatterin^l^. 

For XRS spectra for q _L c, the experimental results 
and first principles calculations have a good qualitative 
agreement on both the general spectral structure and 
peak positions (Fig. 3). However, there is a disagreement 
about the g-dependence of the main edge peak at 199 eV. 
The experimental data shows about 15% increase in this 
peak's intensity when progressing from 2 to 9A , while 
the theory predicted a decreasing trend. The discrep- 
ancy is too large to be accounted for by the randomly 
oriented h-BN crystallites within the basal plane, and 
may be partially mediated by a weak g-dependence in 
the atomic backgrouncPSl. Future studies are needed to 
clarify it. However, this discrepancy does not impact our 
main results. 

We also performed a crude model calculation of this 
exciton 's g-dependence using a hydrogen-like 2pz final 
state with a screened (static dielectric constant e = 7) 
Z + I approximation of the core-hole potential from Eq. 
1 (dashed line in Fig. 4)|£12SI^ resulting in a g-dependence 
much narrower than the experimental data. The failure 
of this calculation suggests a deeply-bound Frenkel-type 
exciton^ , which would provide a broader final momentum 
distribution. 

In conclusion, we have used non-resonant x-ray Raman 
scattering (XRS) to perform exciton spectroscopy on h- 
BN, and find a liQ-type exciton at the pre-edge position 
of the boron K-shell spectra for q || c. This conclusion 
is supported by independent first principles calculations 
of the XRS and of the projected density of states. The 
existence of a pure-type exciton is strongly related to the 
presence of reflection symmetry in the system. This work 
illustrates the role of symmetry-based arguments in the 
interpretation of single-crystal q-dependent XRS studies 
and more broadly illuminates the interplay between lo- 
cal atomic symmetry and final-state effects in core-shell 
spectroscopies. 
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